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Abstract 

The  power  output  and  efficiency  expressions  for  thermoelectric  (semiconductor)  generators  which  is  composed  of  multi-elements  are 
derived  with  considerations  of  heat  transfer  irreversibility  in  the  heat  exchangers  between  the  generator  and  the  heat  reservoirs.  Numerical 
examples  are  provided.  The  effects  of  heat  transfer  and  the  number  of  elements  on  the  performance  are  analyzed.  ©  2002  Editions  scientifiques 
et  medicales  Elsevier  SAS.  All  rights  reserved. 
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1.  Introduction 

Several  authors  have  applied  the  finite  time  thermody¬ 
namics  or  entropy  generation  minimization  [1-3]  to  the 
analysis  of  thermoelectric  generator  [4-10].  They  analyzed 
the  effect  of  finite-rate  heat  transfer  between  the  thermo¬ 
electric  device  and  its  external  heat  reservoirs  on  the  per¬ 
formance  of  the  single-element  thermoelectric  generator.  In 
practice,  a  thermoelectric  generator  is  composed  of  many 
fundamental  thermoelectric  elements.  It  is  a  multi-element 
device.  In  the  present  paper  we  will  investigate  the  charac¬ 
teristics  of  a  multi-element  thermoelectric  generator  with  the 
irreversibility  of  finite-rate  heat  transfer,  Joulean  heat  inside 
the  thermoelectric  device,  and  the  heat  leak  through  the  ther¬ 
moelectric  couple  leg.  The  expressions  for  power  output  and 
efficiency  of  the  multi-element  thermoelectric  generator  are 
derived  algebraically  and  illustrated  numerically. 

2.  Power  and  efficiency  analyses 

Real  thermoelectric  generator  is  composed  of  many 
thermoelectric  elements,  which  is  shown  in  Fig.  1.  Each 
element  is  composed  of  P-type  and  N-type  semiconductor 
legs.  The  thermoelectric  element  is  assumed  to  be  insulated. 


both  electrically  and  thermally,  from  its  surroundings,  except 
at  the  junction-reservoir  contacts.  The  internal  irreversibility 
is  caused  by  Joulean  electrical  resistive  loss  and  heat 
conduction  loss  through  the  semiconductors  between  hot 
and  cold  junctions.  The  Joulean  loss  generates  an  internal 
heat  I2R,  where  R  is  the  total  internal  electrical  resistance 
of  the  semiconductor  couple  and  I  is  the  electrical  current 
generating  from  the  couple.  The  conduction  heat  loss  is 
K  (7\vh  —  7\vl),  where  K  is  the  thermal  conductance  of  the 
semiconductor  couple,  7wh  is  the  hot  junction  temperature, 
and  ?wl  is  the  cold  junction  temperature.  The  external 
irreversibility  is  caused  by  finite  rate  heat  transfer,  i.e.,  the 
temperature  differences  (7h  —  Twh)  and  (7\vl  —  7jJ.  Also 
{K\F\)i  and  {K^F^i  in  Fig.  1  are  the  heat  conductances  in 
the  hot-  and  cold-side  heat  exchangers  for  ith  thermoelectric 
element,  where  K\  and  fj  are  the  heat  transfer  coefficient 
and  the  heat  transfer  surface  area  in  the  hot-side  heat 
exchanger  for  ith  thermoelectric  element,  and  K2  and  Fi 
are  the  heat  transfer  coefficient  and  the  heat  transfer  surface 
area  in  the  cold-side  heat  exchanger  for  ith  thermoelectric 
element.  Therefore,  for  each  element,  we  have  the  rate  (Q'H) 
of  heat  transfer  from  the  heat  source  Tfi  to  the  element 
hot  junction  at  temperature  7wh,  and  the  rate  ( Q[)  of  heat 
transfer  from  the  element  cold  junction  at  temperature  7\vl 
to  the  heat  sink  7p  as,  respectively, 

Gh  =  [«7wh/  -  0.5 12R  +  K(TW H  -  TWl)] 

=  (ATiFiMTh-Twh) 
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Nomenclature 


F\  heat  transfer  surface  area  in  the  hot-side 

heat  exchanger .  m 

F2  heat  transfer  surface  area  in  the  cold-side 

heat  exchanger .  m 

I  electrical  current  generating 

from  the  semiconductor  couple .  A 

/p  optimal  electrical  current  generating  from 

the  semiconductor  couple  at  maximum 

power  point .  A 

Ir)  optimal  electrical  current  generating 

from  the  semiconductor  couple  at  maximum 

efficiency  point .  A 

K  thermal  conductance 

of  the  semiconductor  couple .  W  m_1  K_1 

K\  heat  transfer  coefficient 

in  the  hot-  side  heat  exchanger _  W-  m_  1  •  K- 1 

K2  heat  transfer  coefficient 

in  the  cold-side  heat  exchanger . . .  W  m_1  K_1 
n  number  of  the  thermoelectric 

P  elements  power  output .  W 

Pmax  maximum  power  output .  W 

Pt]  power  output  at  maximum  efficiency  point . .  W 


Qh  heat  transfer  from  the  heat  source 

to  the  hot  junction  at  temperature .  W 

<2l  heat  transfer  from  cold  junction 

to  heat  sink .  W 

R  total  internal  electrical  resistance 

of  the  semiconductor  couple .  Q 

Ri  electrical  resistance  of  the  external  load .  Q 

7h  temperature  of  the  heat  source .  K 

7l  temperature  of  the  heat  sink .  K 

7wh  temperature  of  the  hot  junction .  K 

7wl  temperature  of  the  cold  junction .  K 

Greek  symbols 

o;n  Seebeck  coefficient  of  the  N-type 

semiconductor  leg .  V-  K- 1 

ap  Seebeck  coefficient  of  the  P-type 

semiconductor  leg .  V-  K- 1 

r]  thermal  efficiency 

bmax  maximum  thermal  efficiency 

r]  p  thermal  efficiency  at  maximum  power 

point 


Q'l  =  [«7wl/  +  0.5 12R  +  K(Twli  -  Twl)] 

=  (K2F2)i(TWL-TL)  (2) 

where  a  =  ap  —  aNDap  and  cxn  are  the  Seebeck  coefficients 
of  the  P-  and  N-type  semiconductor  legs. 

For  a  real  thermoelectric  generator  composed  of  n  iden¬ 
tical  thermoelectric  elements,  we  have  its  rate  of  absorb¬ 
ing  heat  (<2h)  and  rate  of  releasing  heat  (Ql)  as,  respec¬ 
tively: 


Th 


Q’4(KlF'f 


£h  =  Ki  F]  (Thc  -TH)=n  Q'n  (3) 

Ql  =  K2F2(Tl  -  7lC)  =  nQ'L  (4) 

where  K\F\  and  K2F2  are  the  heat  conductances  (prod¬ 
uct  of  heat  transfer  coefficient  and  heat  transfer  surface 
area)  of  heat  exchangers  between  the  hot  and  cold  junctions 
of  the  thermoelectric  generator  and  their  respective  reser¬ 
voirs. 

Combining  Eqs.  ( 1 )— (4)  yields: 

7wh  =  ((KlFlTH/n+0.5I2R))(K2F2/n  +  aI) 

+  K(KiF\T}i/n  +  K2F2TL/n  +  I2R ) 
x  ((KiFx/n  +  aI){K2F2/n  -  <xl) 

+  K(KlFl  +  K2F2)/n)~l  (5) 

Twl  =  {{K2F2Th/n  +  0.5I2R))(K1Fl/n  +  aI) 

+  K(KiFiTH/n  +  K2F2TL/n  +  I2R ) 
x  ((KiFx/n  +  aI){K2F2/n  -  al) 

+  K(KlFl  +  K2F2)/n)~l  (6) 

Qh  =  KiFi{l3Ra/2-I2[Tm2  +  K2F2R/(2n)  +  KR] 

+  IK2F2THa/n  +  KK2F2(TH  -  7jJ/n} 
x  ({KxFx/n  +  aI){K2F2/n-Ia) 


Fig.  1 .  Thermoelectric  generator  element. 


+  K(KlFl  +  K2F2)/n)' 


(7) 
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Q  L  =  K2F2{l3Ra/2  +  I2[TLa2  +  KlFiR/(2n)  +  KR ] 

+  IK\F\T\^a/n  +  KK,F,(TH  -  Th)/n } 
x  ((KiFi/n  +  aI)(K2F2/n  —  al) 

+  K(KlFl+K2F2)/n)~1  (8) 

The  total  power  output  ( P )  and  the  efficiency  (/))  of  the 
multi-element  thermoelectric  generator  are  as  follows: 


P  =  Qh  ~  Ql  (9) 

n  =  1  -  Ql/Qk  (10) 

Combining  Eqs.  (7)— (10)  gives: 


P  =  (I3(KiFi  -  K2F2)Ra/t 

-  I2[(KlFlTH  +  K2F2TL)a2  +  K\F\K2F2R/n 

+  (KlFl  +  K2F2)KR] 

+  IKlFlK2F2a(TH-TL)/n) 
x  ((KiFi/n+alKK^/n-al) 

+  K(K1F1  +  K2F2)/n)~1  (11) 

r /  =  (I3(KiFi  -  K2F2)Roi/2 

-  /2[(^iFi7h+  K2F2Th)a2  +  K\F\K2F2R/n 

+  (KlFl  +  K2F2)KR] 

+  IKyFyK2F2a(TH-TL)/n) 
x  (K\F\{l3 Ra/2 

-  I2[THa2  +  K2F2R/(2n )  +  KR ] 

+  IK2F2THa/n 

+  KK2F2(Th  —  7L)/«})-1  (12) 


Eqs.  (11)  and  (12)  are  the  major  results  of  this  paper. 
They  reflect  the  effects  of  heat  transfers  (K\F\  and  K2F2), 
heat  reservoir  temperatures  (7k  and  7l),  internal  heat 
conductance  (K),  internal  electrical  resistance  (R),  Seebeck 
coefficient  (a),  working  electrical  current  (/)  and  number 
of  thermoelectric  elements  ( n )  on  the  power  output  (P)  and 
efficiency  (rf)  of  a  multi-element  thermoelectric  generator. 

If  n  =  1,  Eqs.  (11)  and  (12)  become  the  finite  time 
analysis  results  of  single-element  thermoelectric  generator 
[4,5,9,10]. 

If  K i  F\  —  K2F2  -»■  oo,  7\vh  =  7h  and  7\vl  =  7L, 
Eqs.  (11)  and  (12)  become  theresults  of  conventional  non¬ 
equilibrium  thermodynamic  analysis 


P  =  [«(7h  -  7jJ/  -  I2R]n 
=  a(TH-TL)I-I2R 
11  aTnI  -  0.5I2R  +  K(Th-  Tl) 


(13) 

(14) 


In  this  case,  P  is  dependent  on  n  but  rj  is  independent 
of  n.  If  n  =  1,  Eqs.  (13)  and  (14)  become  the  results  of 
references  [1 1,12],  However,  7wh  =  7h  and  7\vl  =  T\.  need 
infinite  heat  transfer  surface  area,  its  specific  power  output 
P/(F]  +  F2)  is  zero  [7], 


Eqs.  (11)  and  (12)  show  that  there  exist  two  important 
points:  a  maximum  power  point  (Emax)  with  the  correspond¬ 
ing  efficiency  rj p  and  the  optimal  electrical  current  /p,  and 
a  maximum  efficiency  point  (j) max)  with  the  corresponding 
power  output  Ptl  and  the  optimal  electrical  current  It].  From 
the  point  of  view  of  finite  time  thermodynamic  optimization 
(the  compromise  optimization  between  power  output  and  the 
efficiency),  the  design  parameters  optimal  region  of  the  ther¬ 
moelectric  generator  should  be: 

Pn^P^  Pmax  (15) 

W  <  r)  <  max  (16) 

For  the  case  of  infinite  rate  heat  transfer,  7p,  /Jmax ,  r] p,  In, 
Omax  and  Pn  could  be  obtained  analytically  [11].  However, 
for  the  general  case  of  finite  rate  heat  transfer,  they  should 
be  obtained  numerically. 


3.  Numerical  examples 

Taking  7k  =  600  K,  TL  =  300  K,  a  =  2.3  x  10“4  V 
•K-1,  R  =  1.4  x  10“3  ft,  K  =  1.5  x  10“2  W-r'-m"1, 
K\F\  =4  W-K_1  and  K2F2  =  1  W-K”1,  the  power  ver¬ 
sus  current  and  efficiency  versus  current  characteristics  are 
calculated  as  shown  in  Figs.  2  and  3.  These  parameters 
are  of  a  product,  IEC 1-0390 18  thermoelectric  generator, 
of  Shanghai  Jin  Wei  Thermo-electricity  Co.  Ltd.  [15].The 
dashed  lines  represent  the  performance  of  infinite  rate  heat 
transfer  thermoelectric  generator.  Figs.  2  and  3  show  that 
the  heat  transfer  irreversibility  does  affect  the  power  output 
and  the  efficiency  of  the  thermoelectric  generator.  This  ef¬ 
fect  must  be  considered  in  the  performance  analysis.  The 
power  output  with  considering  the  heat  transfer  effect  is 
about  half  of  the  power  output  without  considering  the 
heat  transfer  effect  in  the  case  of  n  =  10.  For  the  fixed 
heat  conductance  (K\Fi  and  K2F2 ),  the  maximum  elec¬ 
trical  current  corresponding  to  zero  power  output  and  zero 
efficiency,  the  optimal  current  corresponding  to  the  maxi¬ 
mum  power  output,  the  optimal  current  corresponding  to  the 
maximum  efficiency,  and  the  efficiency  decrease  with  the 
increase  of  number  (n)  of  thermoelectric  elements.  How¬ 
ever,  there  exists  power  extremal  with  respect  to  number  (n) 
of  thermoelectric  elements,  that  is,  there  exists  an  opti¬ 
mal  n  corresponding  to  maximum  power  output.  There¬ 
fore,  there  exists  a  double  maximum  power  with  the  opti¬ 
mal  current  and  the  optimal  number  of  thermoelectric  ele¬ 
ments. 

In  the  conventional  analysis  without  considering  heat 
transfer  effect,  the  power  output  is  a  monotonic  (linear) 
increasing  function  of  n,  and  the  efficiency  is  indepen¬ 
dent  of  n.  Hence,  the  effect  of  heat  transfer  irreversibil¬ 
ity  increases  when  n  increases  for  the  fixed  heat  conduc¬ 
tance. 

In  order  to  analyze  the  effect  of  heat  conductance  on 
the  performance  of  thermoelectric  generator,  we  calculate 
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1(A) 


Fig.  2.  Power  vs.  current  and  number  of  elements  (S'i-F'i  =  4  W-K  1 , 

f2f'2  =  iw-r1). 


1(A) 

Fig.  4.  Power  vs.  current  and  number  of  elements  (K]Fi  =  1  W-K-1, 
K2F2  =  0.25  W-K-1). 


1(A) 

Fig.  3.  Efficiency  vs.  current  and  number  of  elements  (K\F\  =4  W-K- 1 , 
K2F2=  1  W-K-1). 


1(A) 

Fig.  5.  Efficiency  vs.  current  and  number  of  elements  (K\F\  =  1  W-K-1, 
K2F2  =  0.25  W-K-1). 


the  performance  by  taking  K \  F\  =  1  W-K-1,  K2F2  — 
0.25  W-K-1  and  as  shown  in  Figs.  4  and  5  (the  other 
parameters  are  the  same  as  those  for  Figs.  2  and  3). 
Comparing  Figs.  2  and  4,  and  Figs.  3  and  5,  we  know  that  the 
heat  conductance  affects  the  power  and  efficiency  obviously. 
With  the  decrease  of  heat  exchanger  performance,  the 
maximum  current,  the  optimal  current  corresponding  to  the 
maximum  power  output,  the  optimal  current  corresponding 
to  the  maximum  efficiency  and  the  optimal  number  of 
thermoelectric  elements  corresponding  to  the  maximum 
power  output  decrease. 

The  numerical  results  show  that  the  power  vs.  efficiency 
characteristics  of  real  multi-element  thermoelectric  genera¬ 


tor  is  a  loop-shaped  curve.  This  is  consistent  with  those  of  a 
generalized  irreversible  Carnot  engine  theoretical  model  es¬ 
tablished  by  Chen  et  al.  [13,14], 


4.  Conclusion 

The  results  of  this  paper  show  that  the  heat  transfer  ir¬ 
reversibility  does  affect  the  performance  of  thermoelectric 
generator.  This  effect  must  be  considered  in  the  analysis.  For 
the  real  generator  composed  of  multi-element,  the  number 
of  thermoelectric  elements  affects  the  performance  too.  The 
optimal  current  and  the  optimal  number  of  thermoelectric 
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elements  must  be  choose  in  the  points  of  view  of  the  com¬ 
promise  optimization  between  power  output  and  efficiency 
in  order  to  obtain  the  best  performance. 

The  performance  optimization  of  thermoelectric  genera¬ 
tor  could  be  carried  out  further.  The  optimization  includes 
two  aspects,  i.e.,  the  internal  optimization  and  the  external 
optimization.  The  former  is  to  optimize  the  thermoelectric 
couple  leg  size  of  each  element  to  minimize  (KR).  The  lat¬ 
ter  is  to  optimize  the  distribution  of  heat  transfer  surface 
areas  or  heat  conductance  for  the  two  heat  exchangers  and 
to  optimize  the  external  load  (R\  )  and  internal  resistance 
( R )  matching.  Those  will  be  the  subjects  of  a  future  pa¬ 
per. 
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